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Abstract. Studies of the halo gas in the Milky Way and in nearby spiral galaxies show the 

presence of gas complexes that cannot be reconciled with an internal (galactic fountain) origin 

and are direct evidence of gas accretion. Estimating gas accretion rates from these features 

consistently gives values, which are one order of magnitude lower than what is needed to feed 

the star formation. I show that this problem can be overcome if most of the accretion is in 

fact hidden as it mixes with the galactic fountain material coming from the disk. This model 

OO not only provides an explanation for the missing gas accretion but also reproduces the peculiar 

^^ kinematics of the halo gas in particular the vertical rotation gradient. In this view this gradient 

J— ^ becomes indirect evidence for gas accretion. 
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1. Introduction 

A large amount of fresh gas accretion onto the Milky Way has been advocated for 



decades since the work of Twarog (1980) showed the relative constancy of the Star 

^5 Formation Rate (SFR) in the Galactic disk over the last 10 Gyrs. Roughly constant SFR 

requires constant accretion and this translates today into a need for gas accretion at a 



Q rate of the order ~ 1 Mq yr ^ . 



Since the 1960s, the presence of a large number of neutral hydrogen (Hi) clouds filling 



^ the sky around us and having, on average, negative velocities with respect to the Milky 

i__i Way disk has been recognized (e.g. Hulsbosch 1968). Such clouds, named High Velocity 

Clouds (HVCs), given their largely anomalous radial velocities (Wakker & van Woerden 

1997), were immediately regarded as possible evidence for gas accretion from intergalactic 

space onto our Galaxy (Oort 1970). Alternative interpretations have been proposed, for 

\^ example that the clouds are produced by the cooling of material expelled from the disk 

CO via a "galactic fountain" (Bregman 1980) or that they are a much more distant Local 

CO Group population (e.g. Blitz 1999). 

f-^ From the time of Oort's suggestion until very recently there were two major unknowns 

^D about the HVCs: their distances and their metallicities. It is thanks to a tenacious ob- 

QQ servational campaign carried out especially over the last few years that we now know 

T\ both distances and metallicities for all the major HVCs (Wakker et al. 2007, 2008). The 

^ results leave no doubts: the HVCs are located in the Milky Way halo and have metallic- 

k^ ities of about an order of magnitude lower than the average disk ISM metallicity. These 

'^ properties make it most likely that this is accreting material falling onto the Milky Way 

C^ for the first time. 

The HVCs do not comprise the entire halo population, rather they are an extreme 
population at large heights and having particularly anomalous velocities. In the lower 
halo large amounts of cold gas are observed (e.g. Lockman 1984), some of which go under 
the name of Intermediate Velocity Clouds (IVCs) and are regarded as a likely galactic 
fountain component (e.g. Wakker 2001). More recently, Kalberla fc Pedes (2008) ^ showed 
that up to 10% of the Milky Way HI gas is in fact extra-planar and highly turbulent. 
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2. Halo gas in nearby spiral galaxies 

In the last decade, the study of the halo (extra-planar) gas has been extended to nearby 
galaxies. This is a demanding task, given the low surface brightness of the halo emission, 
and the studies have been restricted to a relatively small number of objects. In Table [l] 
I summarize the results for the best studied galaxies so far. This table includes galaxies 
seen at different inclination angles. For edge-on galaxies the halo gas can be separated 
spatially from the disk gas, whilst for galaxies seen at intermediate inclinations it can be 
separated thanks to its peculiar kinematics. The main kinematic feature of extra-planar 
gas is its decreasing rotation velocity with increasing height from the plane, it is said to 
be "lagging" behind the disk gas. Such a velocity gradient has been estimated for a few 
galaxies (column 9, Table fl]) and it is an important constraint for models of extra-planar 
gas formation (see Section 4) . The presence of this gradient is also the main reason why 
extra-planar gas can be detected in non edge-on galaxies (Fraternali et al. 2002). 

Table 1. Physical properties of extra-planar gas in spiral galaxies 



Galaxy 


Type 


incl 


Vflat 


Mnihaio 


Mnitot 


SFR 


Acer, rate 


Gradient"" 


Ref. 






n 


(km/s) 


(10® Mo) 


(10^ Mo) 


(Mo/yr) 


(Mo/yr) 


(km/s/kpc) 




Milky Way 


Sb 


- 


220 


^4 


4 


1-3 


^0.2^ 


-22 


(1,2,3) 


M31 


Sb 


77 


226 


>0.3 


3 


0.35 


- 


- 


(4,5) 


NGC 253 


Sc 


-75 


-185 


0.8 


2.5 


> 10 


- 


- 


(6) 


M33 


Scd 


55 


110 


>0.1 


1 


0.5 


0.05" 


- 


(7,8) 


NGC 2403 


Scd 


63 


130 


3 


3.2 


1.3 


0.1 


12 


(9) 


NGC 2613 


Sb 


-80 


-300 


4.4^ 


8.7 


5.1 


- 


- 


(10) 


NGC 3044 


Sc 


84 


150 


4 


3 


2.6^ 


- 


- 


(11) 


NGC 4559 


Scd 


67 


120 


5.9 


6.7 


0.6" 


- 


~-10 


(12) 


NGC 5746 


Sb 


86 


310 


- 1 


9.4 


1.2 


0.2^ 


- 


(13,14) 


NGC 5775 


Sb 


86 


200 


- 


9.1 


7.7" 


- 


-8" 


(15) 


NGC 6946 


Scd 


38 


175 


^2.9 


6.7 


2.2 


- 


- 


(17) 


NGC 891 


Sb 


90 


230 


12 


4.1 


3.8 


0.2 


-15 


(18) 


UGC7321 


Sd 


88 


110 


>0.1 


1.1 


-0.01^ 


- 


> -25 


(19) 



" Gradient in rotation velocity with height (from the flat part of the rotation curve); ^ from 
complex C and other clouds with known distances in (2) without correction for the ionised 
fraction; " from the HI mass in (8) without their correction for the ionised fraction; ^ from sum 
of the various ext ra-planar clouds; " calculated from the FIR luminosity using the formula in 



Kewley et al. (2002) , ^ from the counter-rotating cloud (13) using an infall time-scale of 1 x 10^ 
yr; ^ calculated using opt ical lines (16); ^ SFR of only massive stars > 5 M p. References: (1) 
Kalberla k Dedes (2008) » (2) "Wa kker et al. (2007)1 IWakker et al (2008)| (3) 'Levine, Hei lesfc] 



felitz (2008); (4) Thilker et al. (2004 )| (5)|Walterbos fc Braun (1994)1 (6)|B oomsma et al. (2005 



miReakes & Newton (1978)i (8) iGrossi et al . (2008); (9) |Fraternali et al. (2002)[ (10 ) Chaves 
ttlTw in (2001); (11) Lee fc Irwin (1997)| (1 2) Barbi eri et al. (2005)| (13) |Rand fc~B enjamin' 
r2008)|(14)IPedersen et al. (2006); (15) Irwi n (1994)| (16 ) IHeald et al. (2006a)| (17)|Boo msma 
et al. (2008)1 (18) |Oost'erIoo, Fraternali fc Sancisi (2007)| (19) [Mattliews fc Wood (2003) 



Fig.[T]shows the two types of extra-planar gas detections. On the left, the total HI map 
of the edge-on galaxy NGC 891 (blue + contours) is overlaid onto a DSS optical image 
(orange). The HI data (obtained with the Westerbork Synthesis Radio Telescope) show 
a massive and extended H I halo with a mass of 1.2 x 10^ Mq or 30% of the total H I mass 
(Oosterloo, Fraternali & Sancisi 2007). On the right panel of Fig. fl] the position- velocity 
(p-v) plot along the major axis of the galaxy NGC 2403 (inclination 60°) obtained with 
the VLA (Fraternali et al. 2001). A broad component of gas at rotation velocities lower 
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Figure 1. Two methods of detecting extra-planar gas. Left: total HI map (blue + contours) for 
the edge-on galaxy NGC 891 overlaid on the optical image (orange) (data from Oosterloo et al. 
2007). The extra-planar gas is clearly separated on the sky. It surrounds the whole galactic disk 
with a filament that extends up to 20 kpc. Right: posit ion- velocity diagram along the major axis 
of the intermediate inclination galaxy NGC 2403 (from Fraternali et al. 2001). The extra-planar 
gas is kinematically separated from the disk gas. It is seen as a faint component rotating more 
slowly than the disk (white dots = disk rotation curve). 

than the disk gas (also called the "beard") is clearly visible at low emission levels (light 
grey). This beard component is the halo gas in NGC 2403. 

Extra-planar gas is possibly ubiquitous as several nearby galaxies other than those 
reported in Table [T] show indications of its presence. It is also observed in the ionised 
phase. Optical studies of nearby edge-on galaxies show that roughly half of them have 
extended layers of diffuse ionised gas (e.g. Rossa & Dettmar 2003) and with similar 
kinematics as the HI layers (Heald et al. 2006b). 



3. Direct evidence for gas accretion 

Direct evidence of gas accretion is difficult to obtain. The strategy is to look for gas 
components (usually at very anomalous velocities) which are incompatible with an in- 
ternal origin. The large majority of the extra-planar gas studied so far has actually a 
very regular kinematics that follows closely the kinematics of the disk (see for instance 
the p-v diagram for NGC 2403, right panel of Fig. [l]). This points to a tight connection 
between disk and halo components. 

The first features that deserve attention in the search for gas accretion are the HI 
filaments. Among the galaxies in Table[T]at least half show large filamentary H I structures 
in their halos, the most notable cases being NGC 891 and NGC 2403. NGC 891 has a long 
massive filament (Mhi ~ 1.6 x 10^ M©) extending up to about 20 kpc from the plane 
of the disk (Fig. fl]). NGC 2403 also has a filament with a similar HI mass located in 
projection outside the bright optical disk and clearly separated in velocity from the disk 
kinematics (Fig. [2]). They are both very similar to Complex C in our Galaxy. We can 
calculate the energy needed to form these filaments assuming that they come from the 
disk through a galactic fountain. In the case of the NGC 891 filament it turns out that 
this energy should be of the order ~ 1 x 10^^ erg. This would correspond to the explosion 
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Figure 2. The HI filament in NGC2403. Left: optical image overlaid with the HI map (in 
contours) of the 1 x 10^ Mq filament. Right: position- velocity plot along the dashed line in the 
left panel. The filament is 8kpc long and separated from the normal disk kinematics (higher 
contours) by about 90kms~^. Note the small dispersion in velocity of the filament and its 
proximity to the systemic velocity of the galaxy (t'sys = 133kms~^). 

of about 10^ supernovae in a specific region of space and over a time-scale shorter than 
a dynamical time, which is clearly a very unlikely event. 

A second type of feature that has been found in these new deep surveys are clouds 
at very anomalous velocities that end up in the region of counter-rotation. The data 
of NGC 891 show two such clouds with masses of order ~ 10^ Mq and counter-rotating 
velocities of ~ 50 and ~ 90kms~^ (Fig.[3|. These clouds cannot be produced in any kind 
of galactic fountain and they are most likely direct evidence of gas accretion. NGC 2403 
also has gas components at very anomalous velocities called "forbidden gas" (Fig. [ij Fra- 
ternali et al. 2002). Counter-rotating or forbidden clouds are also observed in NGC 4559, 
in NGC 5746 and NGC 6946. 




Figure 3. Counter-rotating clouds in NGC 891. These two channel maps at radial velocities 580 
and 615 kms~^ (vsys = 528kms~^) show gas on the receding side (to the right of the vertical 
dashed line) of the galaxy. Two clouds are visible in isolation on the left side of the galaxy 
(arrows), they are counter-rotating (from Oosterloo et al. 2007. 



If we believe that the above structures have an extragalactic origin we can estimate the 
rate of gas accretion by assuming typical infalling times of few x 10^ — 10^ yr. The resulting 
rates are shown in Table 111 (column 8), they are typically of the order 0.1 Mq yr~ and 
generally 1 order of magnitude lower than the SFRs. These directly observed accretion 
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rates include only Hi, and they should be corrected for helium and possibly ionised gas 
fractions. However, it appears difficult to reconcile them with the rates of star formation 
(column 7, Table IT]). This result is very common for nearby galaxies (Sancisi et al. 2008). 

4. Indirect evidence for gas accretion 

The result that the rate of gas accretion onto galaxies which is directly observed is much 
lower than expected implies that most of the accretion should be somewhat "hidden" . 
I describe here possible indirect evidence of this missing gas accretion, provided by the 
rotation velocity gradient of the extra-planar gas. The steepness of this gradient is not 
reproduced by galactic fountain models (e.g. Fraternali & Binney 2006; Heald et al. 
2006b) as they tend to predict shallower values (a factor half or less). Fig. [4] highlights this 
problem for NGC891. The points are rotation velocities derived at heights z = 3.9 kpc 
and z = 5.2 kpc from the plane (Fraternali et al. 2005). Clearly the fountain clouds in 
the model rotate too fast (have a larger angular momentum) than the extra-planar gas 
in the data. 
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Figure 4. Indirect evidence for gas accretion. Rotation velocities (black points) of the HI 
extra-planar gas in NGC891 at 3.9 kpc (left) and 5.2 kpc (right) from the plane compared to 
the disk rotation curve (dotted line) and predictions from two models: a pure galactic fountain 
model (red line above) and a model where the fountain clouds sweep up and accrete ambient gas 
during their passage through the halo (blue line below). The accretion rate required to produce 
this fit is ^ 3M0yr~^, very similar to the star formation rate of NGC891 (from Fraternali & 
Binney 2008). 



How can the fountain clouds loose part of their angular momentum? Fraternali fc 
Binney (2008) | consider the possibility that fountain clouds sweep up ambient gas as 
they travel through the halo. In this scheme ambient gas condenses onto the fountain 
clouds, these latter grow along their path through the halo and eventually fall down into 
the disk (see Fig. |5|. If the ambient gas has relatively low angular momentum about 
the z-axis then this condensation produces a reduction in the angular momentum of the 
fountain gas. The only free parameter of the model is the accretion rate, which is tuned 
to reproduce the rotation curves of the extra-planar gas. Remarkably, the required gas 
accretion rate turns out to be very similar to the SFR. For NGC 891 we found a best- 
fit accretion rate of about 3 Mq yr~^ (see the blue curves in Fig. |4| and for NGC 2403: 
0.8 Mq yr~^. In NGC 2403, this model is also able to reproduce the observed radial inflow 
of the halo gas (see Fraternali & Binney 2008). 

One implication of the above fountain+accretion model is that it predicts that most 
of the extra-planar gas is produced by the galactic fountain and only a small fraction 
(about 10%) is extragalactic. This is in agreement with the metallicity of the IVCs and 
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Figure 5. Schematic of th e fountain gas sweeping up ambient gas. The model presented in 
Fraternali & Binney (2008) suggests that superbubbles blow gas into the halo which then sweeps 
up part of the ambient gas during its passage before falling back to the disk. This mechanism 
allows fountain clouds to loose part of their angular momentum as required by the data (see 
Fig. p| and produces a net gas accretion at a rate similar to the SFR. 

the clear links between anomalous velocity clouds and star forming regions (e.g. Boomsma 
et al. 2008). A second implication is that it does not require that the accreting gas is in 
any particular phase but only that its angular momentum about the z-axis is less than 
about half the angular momentum of the disk material. Finally, this model predicts an 
accretion rate of the order of the SFR and in general, proportional to the supernova 
rate. Interestingly, this appears to be a general requirement for galaxies throughout the 
Hubble time as it reconciles the observed cosmic star formation history with the gas mass 
in galaxies at low and high redshifts (Hopkins, McClure- Griffiths & Gaensler 2008). 
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